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Abstract 

A random vortex world-surface model for the infrared sector of SU (4) Yang- 
Mills theory is constructed, focusing on the confinement properties and the behav¬ 
ior at the deconfinement phase transition. Although the corresponding data from 
lattice Yang-Mills theory can be reproduced, the model requires a more complex 
action and considerably more tuning than the SU{2) and SU{ 3) cases studied 
previously. Its predictive capabilities are accordingly reduced. This behavior has 
a definite physical origin, which is elucidated in detail in the present work. As the 
number of colors is raised in Yang-Mills theory, the corresponding infrared effec¬ 
tive vortex description cannot indefinitely continue to rely on dynamics determined 
purely by vortex world-surface characteristics; additional color structures present 
on the vortices begin to play a role. As evidenced by the modeling effort reported 
here, definite signatures of this behavior appear in the case of four colors. 
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1 Introduction 


The vortex picture of the strong interaction vacuum was originally conceived [1-5] as 
an explanation for confinement. More recent work [6-23] has established it as a com¬ 
prehensive, and in many respects quantitative infrared effective model of the vacuum, 
generating not only confinement, but also the spontaneous breaking of chiral symmetry 
and the axial Ua{ 1) anomaly. The evidence for the vortex picture rests on two pillars. 
On the one hand, methods have been developed which permit the identification and 
isolation of vortex structures in the gauge configurations of lattice QCD [6,7,19]. Based 
on these techniques, the phenomenology induced by the vortices present in lattice gauge 
configurations can be studied. Highlights of the results obtained in this manner include 
vortex dominance of the string tension [6,7,19], the absence of nonperturbative effects 
(string tension, chiral condensate, topological susceptibility) when vortices are removed 
from lattice gauge configurations [18,19], and a natural explanation of the deconfining 
phase transition as a vortex percolation transition [16]. 

On the other hand, an infrared effective random vortex world-surface model has 
been developed [21-27] which demonstrates that the main nonperturbative features 
of the strong interaction vacuum can be described 2 on the basis of a weakly coupled 
vortex dynamics. This further buttresses the notion that vortices constitute the relevant 
infrared gluonic degrees of freedom in that vacuum. Through its reduced set of degrees of 
freedom and the correspondingly simplified dynamics, the random vortex world-surface 
model has allowed the investigation of a wider range of scenaria than has been accessible 
using the vortices extracted from lattice gauge configurations discussed further above. In 
the random vortex world-surface model, not only has the SU (2) case been investigated 
comprehensively [21-23], but also the confining properties of the 577 (3) case have been 
studied in detail [24-26], and the present work extends the model to 577(4) color. 


2 Motivation 

2.1 Vortex color structure 

To clarify the motivation for the present investigation, it is useful to begin by reviewing 
the different color structures which can occur on center vortices, specifically in various 
Abelian formulations. As will be discussed below, vortices are not only characterized 
by their location in space-time, but possess considerable additional structure in color 
space. In general, it cannot be excluded that this structure may significantly influence 
vortex dynamics; besides geometrical world-surface characteristics, the vortex effective 

2 In the SU{2) and SU( 3) models studied to date, adjustment of a single dimensionless coupling con¬ 
stant suffices to arrive at a phenomenologically viable model, a number of nonperturbative observables 
being quantitatively predicted as a result. The SU( 4) case investigated here turns out to behave less 
favorably in this respect, and the reasons behind this represent a principal focus of the present work. 
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action potentially may also include a dependence on the aforementioned color structure. 
Investigation of the possibility or even necessity of such a dependence is a principal focus 
of the present investigation. Thus, a review of vortex color structure is in order. 

Center vortices are closed tubes of quantized chromomagnetic flux in three spa¬ 
tial dimensions. Accordingly, they are described by (thickened) world-surfaces in four¬ 
dimensional (Euclidean) space-time. The quantization is defined by the property that 
a Wilson loop encircling a vortex (more precisely, linked to the latter) acquires a factor 
corresponding to a center element of the gauge group. In SU(N) gauge theory, there are 
thus N — 1 possible vortex fluxes (the trivial unit element corresponds to no flux being 
present). In particular, this entails that vortices can branch for N > 3, as described in 
greater detail below. 

This description of vortices in terms of their space-time location and their influ¬ 
ence on Wilson loops is in some respects incomplete. While it is adequate for the 
discussion of the confinement properties encoded in Wilson loops, to fully describe the 
topological properties of vortices and the related chiral symmetry breaking phenom¬ 
ena, it is necessary to take into account the direction of the vortex field strength in 
color space [12-14,22,23,28,29]. While this direction can be rotated by gauge trans¬ 
formations, certain properties of it carry gauge-independent significance. For instance, 
generic vortex world-surfaces are nonorientable; this implies that the field strength can¬ 
not be globally proportional to a constant color vector. Instead, the color direction of 
the field strength must vary, in any gauge. While the detailed local space-time form of 
the variation forced by vortex nonorientability is gauge-dependent, it contains invariant 
characteristics which manifest themselves, e.g., in the topological charge. 

Often, it is convenient to construct vortex configurations in an Abelian gauge. In this 
case, variations of the field strength color vector are compressed into lines on the vortex 
world-surfaces, leaving the color vector constant elsewhere. At these lines, the held 
strength is discontinuous in a way which corresponds to the presence of a source or sink of 
magnetic flux, i.e., an Abelian magnetic monopole. Thus, due to the nonorientability of 
vortex world-surfaces, monopoles represent an intrinsic feature of vortex configurations 
in Abelian gauges. 

From this more detailed perspective, in particular also the branching points (in 
three-dimensional space) or lines (in four-dimensional space-time) of vortices mentioned 
further above can be characterized more thoroughly. Vortex branchings in general are 
associated with nontrivial rearrangements in the vortex color structure; they cannot be 
realized while keeping all vortex held strengths in the branching region pointing into 
the same hxed color direction, as will be discussed below. 

To facilitate this discussion, (Abelian) vortex color direction will be specihed in the 
following by giving the (diagonal) chromomagnetic hux matrix Q which results when 
evaluating a closed line integral encircling the vortex, 



This eliminates the detailed space-time information on the vortex field strength which 
is non-essential to the following discussion. In terms of Q, the Wilson loop along the 
same path is 

W = ^ Tr exp(iQ) . (2) 

Note that every vortex world-surface can be associated with two possible orientations, 
which yield opposite signs for Q. Conversely, of course, inverting the integration direc¬ 
tion in the line integral defining Q inverts the sign of the latter. In order to preserve 
rotational symmetry in the following, if a given description allows for vortices associated 
with the color matrix Q, it will also include those associated with the color matrix —Q. 

2.2 Minimal Abelian construction 

The color structure of vortex branchings can now be viewed in a variety of ways. One 
possible picture results if one chooses a minimal set of definite (Cartan) color directions 
to describe the N — 1 types of vortex flux permitted by the gauge group. Under this very 
restrictive condition, vortex branching is tied to the presence of magnetic monopoles. 
Consider the SU (3) gauge group as an example. The two nontrivial center elements are 

exp(27ri/3) , exp(— 2ni/3) . (3) 

These are the phases which can arise in a Wilson loop due to linking with vortices. 
Since these two phases are complex conjugates of one another, they can actually be 
generated by the two possible orientations of one single type of vortex. Consequently, if 
one insists on using a minimal set of color directions, only one color matrix Q , along with 
its negative —Q, is necessary to construct all possible vortex configurations. Without 
loss of generality, let 

Q = y diag (1,1, -2) . (4) 

Now, consider a vortex branching. The only way to divide an incoming vortex flux Q 
into two outgoing fluxes is by letting the latter each be associated with the flux —Q. 
Otherwise, evaluating Wilson loops encircling the incoming vortex or both the outgo¬ 
ing vortices, respectively, would yield different results, violating the Bianchi constraint 
(continuity of flux modulo 27 t). However, even in the allowed branching configuration, 
the incoming flux is not simply the sum of the outgoing fluxes; the flux matrix is dis¬ 
continuous. There is a sink of flux 3 Q at the branching (which is compatible with the 
Bianchi constraint); this is an Abelian magnetic monopole 3 . 

3 Note furthermore that, in this particular case, branchings are the only locations where monopoles 
can be present. A switch from a flux Q to a flux —Q along a single vortex does not correspond to 
a monopole allowed by the Bianchi constraint. This is different from the 5/7(2) case, where a switch 
from Qsu( 2 ) = 7rdiag(1,-1) to —Qsu{ 2 ) is possible. In the minimal description of 5/7(3) vortices, 
therefore, the vortex branching and the Abelian monopole concepts can be used synonymously. Note 
that this strict identification does not persist for 5/7(4); in that case, monopoles can be present away 
from branchings even in the minimal description. 
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As a consequence, one can also take an alternative stance. Instead of discussing 
branchings for N > 3 in particular, and noting that monopoles are present at branchings, 
one can take the point of view that vortex configurations in general contain monopoles, 
which are sources or sinks of up to N vortex fluxes. Monopoles are connected by vortices, 
and between monopoles, vortex fluxes are oriented. Note that, in this language, there is 
no conceptual distinction between the cases N = 2 and N > 3. Both cases are described 
in the same language. In the former case, only two vortices emanate from any monopole, 
whereas for higher N, there can be two or more vortices emanating from monopoles, 
cf. also Fig. ffl below. This subsumes the statement that vortices branch for N > 3, 
while maintaining the same conceptual framework for all N. 

2.3 Nexus construction 

Alternative ways to describe the color structure of vortex branchings result if one allows 
for a non-minimal set of color flux matrices. In particular, if one uses a sufficiently 
large set of different flux matrices Q along with their negatives —Q, one can completely 
disassociate Abelian monopoles from branchings. Such a description of branchings was 
introduced in [12-14] along with the term “nexus” to refer to these objects 4 . Staying 
with the SU( 3) example, the relevant flux matrices are 

y diag (1,1, -2) , y diag (1, -2,1) , y diag (-2,1,1) (5) 

(along with their negatives). By letting each of the three vortices meeting at a branching 
be associated with a different flux matrix, monopoles can be completely avoided at 
branchings 5 . Of course, the global structure of the vortex world-surfaces will in general 
force monopoles to be present elsewhere (which is consistent with the nexus description, 
in contradistinction to the minimal description discussed above). 

On the other hand, note that the nexus language distinguishes the N = 2 and the 
N > 3 cases on a qualitative level. The former case does not allow for (quasi-Abelian) 
nexi, i.e., branchings, whereas all other cases do [12], 

2.4 Dynamical issues 

Whether one views vortex color structure in terms of Abelian monopoles alone (as dis¬ 
cussed at the end of section IT^ll or chooses to differentiate between nexi (i.e., branchings) 

4 More precisely, these are referred to as “quasi-Abelian” nexi in [12-14]. On the other hand, the 
“fully non-Abelian” nexi of [12-14] correspond, after a (singular) gauge transformation, to the “Abelian 
monopoles” of the present work. 

5 This property is likely to single out the quasi-Abelian nexus construction of branchings as the one 
most suited in practice for the purpose of evaluating the topological charge of general SU(3) vortex 
world-surface configurations; this is currently under investigation. 


4 



and monopoles, these objects embody a nontrivial color structure which suggests an al¬ 
ternative characterization of the Yang-Mills vortex vacuum. Instead of viewing vortex 
world-surfaces as the only dynamical degrees of freedom, which happen to force the 
presence of certain color structures through their geometrical properties (nonorientabil- 
ity, branching), one can alternatively view monopoles (and nexi) as bona fide degrees 
of freedom, which, at least to a certain extent, determine the form of the vortex fluxes 
emanating from them and connecting them. Which of these pictures is the more ap¬ 
propriate one is a dynamical question; phrased more formally, is the vortex effective 
action dominated by terms involving geometrical vortex world-surface characteristics or 
by terms involving monopole (and nexus) characteristics? 

In the 577(2) and SU( 3) cases studied to date, effective vortex actions based purely on 
world-surface characteristics appear entirely adequate to reproduce the infrared sectors 
of the respective Yang-Mills theories on a quantitative level. Specifically, vortices are 
controlled by world-surface curvature in the SU (2) and S77(3) models. However, it would 
be rather surprising if this picture persisted indefinitely as the number of colors N is 
increased; on the contrary, a shift towards dynamics influenced by monopole (or nexus) 
characteristics is likely, as has been argued previously [30]. This is expected because, as 
the number of colors rises, center flux can be quantized in ever smaller units, while the 
monopoles always constitute sources or sinks of flux 27r, or multiples thereof (in each 
diagonal color component), cf. Fig. [U Thus, it seems plausible that monopoles become 
the dominant carriers of inertia in vortex configurations, and that the vortex world- 
surfaces then adjust to the monopole positions rather than the latter being determined 
by the world-surface dynamics. This should manifest itself in an effective vortex action 
in which the monopoles (or nexi) attain their own dynamical significance. 



Figure 1: Vortex configurations contain monopoles, cf. section POl For N colors, up to 
N vortex fluxes can emanate from any given monopole. Examples for N = 2 (left) and 
N = 8 (right) are depicted. In order to display a generic case, for N — 8 also a vortex 
associated with center phase e*® = e l7r / 2 (the double line) has been included (all the 
other vortices carry center phase e*® = e m / 4 ). As N rises, center flux can be quantized 
in ever smaller units, while monopoles remain sources or sinks of flux 27T, or multiples 
thereof (in each diagonal color component). 

To search for indications of this shift in the dynamical characteristics, and thus probe 
the limits of applicability of pure random world-surface dynamics for center vortices, 
was the principal motivation for the present 57/(4) work. Indeed, of the effective vortex 
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actions investigated below, only ones including an explicit dependence on branching 
properties prove to be phenomenologically viable. Before proceeding to describe the 
investigation in detail, it should be emphasized that the possibility of such a shift in 
the dynamical characteristics does not imply that the vortex picture as a whole may 
become inappropriate at larger N and that, e.g., an Abelian monopole Coulomb gas may 
become an appropriate description. On the contrary, for any finite N, chromomagnetic 
flux must be constricted into center vortices in order to, e.g., correctly account for 
the lack of confinement of adjoint color sources [8]. The aforementioned shift in the 
characterization of the vacuum merely concerns the specific form of the effective vortex 
action; vortices nonetheless continue to represent relevant infrared degrees of freedom. 
All that is implied is that, at larger N , the Abelian monopoles which intrinsically reside 
on generic non-orientable vortex world-surfaces (in Abelian gauges) may cease to be 
completely dependent objects with no dynamical significance of their own; instead, 
their characteristics may begin to play a role in determining overall vortex configuration 
dynamics 6 . 

3 SU(4) vortex model 

3.1 Degrees of freedom 

Vortex flux quantization in the SU (4) vortex model is determined by the three nontrivial 
center elements i,—i and —1 of the SU( 4) gauge group. The former two elements 
are complex conjugates of one another and thus can be generated by the two possible 
orientations of one type of vortex. Therefore, the SU( 4) model contains in all two 
physically different types of vortices; the ones generating a phase —1 and the ones 
generating the phases ±i depending on their orientation in space-time. 

The SU( 4) random vortex world-surface model is constructed in complete analogy 
to the SU( 2) and SU( 3) cases studied previously. The reader is referred to [21,24] 
for details regarding the physical interpretation of the construction. The vortex world- 
surfaces are modeled by composing them of elementary squares on a hypercubic lattice. 
Each elementary square in the lattice is associated with a value 

q^w(x) E {-1,0,1,2} , q vli (x) = ~q^(x) (6) 

where the square extends from the point x into the positive /j and u directions. The right- 
hand relation is simply a reminder of the behavior of flux under space-time inversions, 
i.e., when the orientation of the vortex surface is reversed. In practice, recording only 
( hw( x ) for /i < v is sufficient. The value q IJV {x) = 0 corresponds to no vortex flux being 
present on the square in question; nonzero values of q^ix) could, e.g., label the first 

6 Merely at infinite N, where the center of the SU (N —> oo) gauge group becomes Z(N —> oo) ~ U (1), 
a pure monopole description may be feasible. There, center quantization of flux ceases to constitute a 
constraint, since arbitrarily small units of center flux are possible. 
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element of the flux matrix Q^{x) associated with the elementary square in a minimal 
Abelian construction, 

Q»u(x) e { ^ cliag (—i, — i, — i, 3) , |diag(l,l, 1,-3) , | diag (2, 2,-2,-2) J (7) 

(since, in the following, only Wilson loops and action densities will be studied, it is not 
necessary to distinguish between the last flux matrix in m and the one corresponding to 
inverse orientation; for a study, e.g., of topological properties, this would be necessary). 

An important point to be noted is that the lattice spacing in this approach is a 
fixed physical quantity implementing the notion that vortices possess a finite transverse 
thickness and must be a minimal distance apart to be distinguished from one another'. 
This is discussed in detail in [21,24], While in the SU(2) and SU(3 ) models, there is 
only one type of center vortex (up to orientation), in the present SU( 4) case, there are 
two physically distinct types of vortices. In general, these can have different thicknesses. 
However, there is no straightforward way of implementing variable thicknesses in the 
present hypercubic lattice formulation of the model; all vortices will be described using 
a single lattice. This model restriction should be kept in mind. 

Ensembles of random vortex world-surface configurations on lattices are generated 
using Monte Carlo methods. The elementary update of a given vortex surface configu¬ 
ration can be effected such that the Bianchi constraint is respected at every step. This 
is achieved by updating all six faces of an elementary three-dimensional cube in the 
lattice at once, in a way which corresponds to adding a vortex flux of the shape of the 
cube surface to the flux previously present. Formally, if the elementary cube in ques¬ 
tion extends from the lattice site x into the positive /i, v and A directions, then update 
simultaneously 

%v(x) —> {<hiv{x) + w) mod 4, q^(x + e x ) -»■ {q^{x + e x ) - w ) mod 4 

qv\(x) —> (q uX (x) + w) mod4, q vX (x + e M ) —> ( q vX (x + e M ) - w) mod4 (8) 

q\^(x) -> (q x ^(x) + w) mod 4 , q Xfl (x + e v ) —>• (q XfJ ,(x + e v ) - w) mod 4 

where the value w G { — 1,1, 2} characterizing the superimposed flux in practice is chosen 
at random with equal probability, and the modulo operation is to be carried out such that 
the result again satisfies q G {—1, 0,1, 2}. Since the update effects a linear superposition 
of two fluxes which satisfy the Bianchi constraint, the updated configuration again 
satisfies that constraint (the modulo operation merely generates shifts by 2tt in the flux 
matrix, which are allowed by the Bianchi constraint; in this way, monopoles appear in 
the vortex configuration as the updates proceed). 

7 In the SU( 2) and SU(3) models, the lattice spacing turns out to be 0.39 fm if one fixes the scale 
by setting the zero-temperature string tension to a(T = 0) = (440 MeV) 2 . 
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3.2 Action 


Finally, the action weighting vortex world-surface configurations must be specified. The 
most general form of the action used is the sum 

S[Q] Squtv [q] + ^area [?] T 'S'branch[?] (9) 

with the individual terms specified as follows. 

The immediate generalization of the dynamics studied in the 577(2) and SU( 3) 
models is a world-surface curvature action, adapted to the present case, in which there 
are two physically distinct types of vortices, 

5'curv[?] = Xm q^\( x )\) + c \\q^( x )q^( x - e\)\) ( 10 ) 

X fl L v<\ 

+C(\q lll/ (x - e v ) q^\{x)\) + C(\q^(x - e v ) q^xix - e A )|)) 

with 

<7(0) = 0 C{ 1) = c n (7(2) = c 12 C( 4) = c 22 . (11) 

Note that this action term (as well as all others below) treats q = ±1 fluxes symmet¬ 
rically, as it should, since these two cases correspond to the two possible orientations 
of the same physical vortex type. As is clear from the expression given, for each link 
extending from the site x into the positive /i direction, all pairs of elementary squares 
attached to that link are examined, excluding those in which the two squares lie in 
the same plane. Depending on the flux associated with the squares, an action incre¬ 
ment is incurred. Thus, world-surfaces are penalized for going around a corner, i.e., for 
curvature. In the following, the couplings c l3 will be generated from two independent 
coefficients ci,c 2 (corresponding to the two physically distinct types of vortices) as 

(‘■/jj CiCj . (12) 

Note that giving c 12 special treatment, i.e., generalizing to c 12 ^ CiC 2 , to an extent would 
amount to a special treatment of branchings (along with related features such as different 
vortex types intersecting along a whole line), since branchings are locations where a 
q — 2 vortex splits into two q — ±1 vortices. However, on the other hand, an additional 
action term weighting branchings explicitly will be introduced separately below. For 
this reason, the option c\ 2 ^ CiC 2 was not pursued further in this investigation (even 
though, strictly speaking, it is physically distinct from the aforementioned branching 
action). Instead, the explicit branching action term introduced below will be studied 
extensively. 

A further type of action investigated in the 577(2) and SU( 3) models is an action 
simply weighting world-surface area, 

5area[g] = A ( I ^ 0*0 I) • 

X fl<U 
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( 13 ) 



In the SU( 2) and SU( 3) models, it was found that this action term could be traded 
off against the curvature term, and it was therefore ultimately discarded. Essentially, 
in the random surface ensemble, total surface area is strongly correlated with total 
curvature content, such that, for all practical purposes, the area term can be replaced 
by a strengthening of the curvature term for a wide range of parameters [21,24], In the 
present work, the effects of this term are briefly explored, cf. section ET21 below: inclusion 
of this term however does not appear to increase the phenomenological flexibility of the 
model. 

Finally, a new term introduced in the present investigation is an explicit weighting 
of branchings, 

‘-’branch[?] = Y Y B ( Y E ^^ X )) + E (%v( X ~ e v)) ] ( 14 ) 

where 

E( 0) = 0 , E — 1 else (15) 

B( 3) = B( 5) = — b , B = 0 else (16) 

This action term precisely encourages vortex branchings: For each link extending from 
the site x into the positive /i direction, the number of attached elementary squares occu¬ 
pied by vortex flux is counted and, if that number equals 3 or 5 (which happens precisely 
when branching takes place), the action is decremented by the branching coefficient b. 
Further below, it will be shown why the introduction of this type of term is necessary 
to arrive at a phenomenologically viable model, corroborating the general arguments 
of section [21 In view of the different vortex formulations discussed in section [21 one 
can interpret Stanch as favoring the presence of nexi, or alternatively as favoring the 
presence of (certain types of) Abelian monopoles. Whichever picture one may choose, 
the introduction of this term represents a departure from pure vortex world-surface dy¬ 
namics, according additional color structures present on the vortex world-surfaces their 
own dynamical significance. 

4 SU(4) Yang-Mills lattice data 

The SU (4) lattice data which will be matched in the random vortex world-surface model 
are drawn from [31,32], After the bulk of the numerical work in the present investigation 
was completed, updated lattice data became available [33]. The author did not attempt 
to adjust the vortex model to the revised data, since the quantitative adjustments this 
would entail do not alter the conclusions drawn. 

On the other hand, [33] also contains lattice data on the spatial string tension, for 
which no measurements were available previously. These data will be compared to 
corresponding predictions within the vortex model constructed here in section 15.31 
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Specifically, the SU (4) Yang-Mills characteristics which will be matched in the vortex 
model are the following: 

4.1 String tension and deconfinement temperature 

The ratio of the deconfinement temperature T c to the zero-temperature quark string 
tension oq is taken to be [31] 

T 

—= = 0.62 . (17) 

VoT 

The quantity Tjy/a was also used in the 517(2) and 517(3) models to £x the single 
curvature coefficient c; since that curvature coefficient represented the only (physically 
significant) adjustable parameter, no further lattice data were needed to fully define 
those models. 

4.2 String tension ratios 

In correspondence to the two different types of center vortices present in the 517(4) 
theory, there are also two distinct string tensions, the quark string tension cq and the 
diquark string tension < 72 . Data on both is needed to fix the two curvature coefficients 
Ci, C 2 of the 517(4) vortex model. Thus, as one generalizes the vortex model to a higher 
number of colors, it does not generally predict ratios between different string tensions, 
unless one makes further prior assumptions about relationships between the properties 
of different types of vortices; for such considerations at large N, cf. [30]. Rather, the 
vortex model always provides sufficient freedom to match all relevant string tension 
ratios, since there are correspondingly many physically distinct types of vortices with 
their separate dynamical characteristics such as curvature coefficients, etc. 

The ratio of the diquark string tension to the quark string tension at zero temperature 
is taken to be [32] 

— = 1.36 . (18) 

Ol 

4.3 First order character of the deconfinement transition 

The first-order character of the deconfinement phase transition becomes more pro¬ 
nounced as one raises the number of colors. In units of T c 4 , the latent heat of 517(4) 
Yang-Mills theory is larger than the one of 517(3) Yang-Mills theory by a factor of 
2.25 [31]. As will be described below, the 517(4) vortex model with a pure curvature 
action m does not reproduce this property, even on a qualitative level. This is the 
reason for the introduction of the new branching action (HI; a more detailed motivation 
for that particular choice follows further below. Introducing this additional action term 
introduces the new branching coefficient 6; to fix this coefficient, an additional piece of 
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lattice data is needed. The aforementioned relation between latent heats in the SU (3) 
and SU(4:) cases can be restated as 


rp 4 
0i 


SU{ 4) 


rn 4 

= 2- As^- 
a 4 


SU( 3) 


(19) 


using (ED as well as the relation T c j a/ct = 0.63 for SU( 3) used in [24], where As denotes 
the difference in action density (per space-time volume) between the two coexisting 
phases at the deconfinement temperature. This modified choice of units was adopted 
because it has more moderate scaling properties as one investigates the vortex model 
on lattices of different (Euclidean) temporal extent N t , cf. below. 

In order to be able to model the property ED, it is necessary to briefly revisit 
the SU( 3) vortex model [24]. The case of lattices with an extension of N t = 2 lattice 
spacings in the Euclidean time direction is discussed in detail in [24]; the deconfinement 
transition is found at the value c = 0.2359 for the curvature coefficient. Reading off 
the discontinuity in the action density at the phase transition from the corresponding 
action density distribution 8 , one has (with a denoting the lattice spacing) 


As = 0.45 T c 4 ( T c a = 1/2, c = 0.2359, aa 2 = 0.695) (20) 


where the value for the zero-temperature string tension a measured separately for c = 
0.2359 has been included for later reference. The physical point of the SU(3) model is 
not at c = 0.2359, but nearby at c = 0.21. This can be corrected for by also considering 
N t — 1 and interpolating in c. For N t — 1, the deconfinement transition is found for 
c = 0.0872. Producing an action density distribution analogous to the one discussed 
above for N t — 2 permits one to read off 

As = 0.022 T c 4 (! T c a = 1, c = 0.0872, era 2 = 1.355) . (21) 


The zero-temperature string tension a was again measured separately for this value of 
c. From ED and ED, the strong variation of As as a function of c in the units given 
is evident. As a consequence, these units are not well suited for interpolation in c, and 
it is indeed advantageous to convert to 



e=0.2359 


0.0075 , 



c=0.0872 


0.0065 . 


( 22 ) 


8 Note that the distributions in [24] are plotted as a function of the action per link divided by the 
curvature coefficient c; therefore, to obtain As, one needs to multiply by 4c/a 4 , where a denotes the 
lattice spacing. Note also that the rough measure for the action density discontinuity used here and 
throughout this work is simply the distance between the two maxima in the action density distribution 
plot (instead of, e.g., modeling the latter as a superposition of two distributions). Since only ratios 
between discontinuities are of interest here, the effects of this slight underestimate cancel to a sufficient 
degree for the present purposes. 
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Finally, interpolating linearly to the physical point c = 0.21, one has 
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and, inserting into PD- the condition to be satished by the SU(4:) model is therefore 
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Figure 2: Distribution of the action density s (per space-time volume) at the deconhning 
phase transition, for a pure vortex world-surface curvature action with c\ = Co = 0.4722. 
The measurement was taken on a 14 3 x 2 lattice; a denotes the lattice spacing. 


5 Exploration of coupling parameter space 

5.1 Pure world-surface curvature action 

As a starting point, consider the SU(A ) random vortex world-surface model with a pure 
curvature action (HDD,©, restricted to equal curvature coefficients, C\ = c 2 . This is 
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not yet very realistic as far as the condition m is concerned; for Ci = c 2 , one obtains 
o\ = uo, i.e., er 2 /<7i = 1. On the other hand, one does observe a deconfinement transition 
which is quite strongly first order, cf. Fig. [21 for N t = 2, taken at C\ = c 2 = 0.4722. 
Similar observations can be made at other N t , with corresponding c\ = c 2 ; the case 
N t = 2 is the one which comes closest to satisfying the remaining condition (1171) . Note 
that the physical set of parameters in the random vortex world-surface model is generally 
such that the inverse deconfinement temperature is not an integer multiple of the lattice 
spacing; properties of the deconfinement transition at the physical point then cannot 
be measured directly, but are inferred by interpolation of data obtained at different N t , 
cf. the treatment of the 577(3) case above and cf. also [21,24], 
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Figure 3: Distribution of the action density s at the deconhning phase transition, for 
a pure vortex world-surface curvature action with c± = 0.4537 and c 2 = 0.5037. The 
measurement was taken on a 26 3 x 2 lattice. 

As a next step, in order to approach the condition m , it is necessary to increase 
c 2 over c\. For example, the set of curvature coefficients Ci = 0.4537, c 2 = 0.5037 
yields <t 2 /cti = 1.17 while simultaneously realizing the deconfinement transition on a 
N t — 2 lattice, cf. Fig. 01 Examining the trends, one notices an unexpected behavior: 
As one attempts to modify the parameters to approach the condition m , the first-order 
character of the phase transition is weakened. In fact, when m is finally satisfied in 
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this way, the action discontinuity at the phase transition has either become undetectably 
small or has disappeared entirely, rendering the deconfinement transition second order. 
To be specific, the author used string tension measurements on N t = 2 lattices in 
addition to the zero-temperature measurements to locate the set of parameters C\ = 
0.428, c 2 = 0.540 which yields CHI) as well as realizing the deconfinement transition at 
N t = 2; then, a careful scan of C\ in steps of 0.0001 around Ci = 0.428 was performed 
(with c 2 = Ci + 0.112 in each case), generating the action density distributions on large 
(30 3 x 2) lattices. No evidence of first order behavior was found. 

This qualitative behavior can be confirmed for other N t . For iV* = 1, the deconfine- 
ment transition occurs at c\ = 0.253, c 2 = 0.353 (with (fTHj) also satisfied). Likewise, for 
N t = 3, the transition occurs at c\ = 0.538, c 2 = 0.571. In both cases, a careful scan in 
Ci, as above, with c 2 — ci fixed, reveals no first-order behavior. By extension, also the 
phase transition at the interpolated set of parameters realizing the remaining condition 
(im> is not discernibly first order. 

Thus, it must be concluded that the conditions m and (01) are incompatible as 
long as one restricts oneself to a pure world-surface curvature action. One can realize 
either the one or the other condition, but one cannot model both simultaneously. 

Before seeking to remedy this by generalizing the action, it is nevertheless useful to 
record here, for later reference, the simplified model which results if one disregards the 
first-order character of the deconfinement phase transition, i.e., if one disregards (1241) . 
and uses a pure world-surface curvature action. The two curvature coefficients Ci and c 2 
are then fixed using only the two conditions CHD and m as follows. Tabled summarizes 
again the pairs of curvature coefficients at which the deconfinement transition occurs 
for N t = 1, 2, 3 while (fTHl) is simultaneously satisfied. 



Cl 

c 2 

T c /y/o 7 

N t = 1 

0.253 

0.353 

0.85 

N t = 2 

0.428 

0.540 

0.60 

£ 

CO 

0.538 

0.571 

0.50 


Table 1: Sets of coupling constants realizing the deconfinement temperature as well as 
satisfying m ■ Condition m is then satisfied by interpolating the parameters using 
the data in the final column. 

The table furthermore records the ratio Tj\[(7\ corresponding to each parameter 
set. To obtain the physical point, one constructs the quadratic interpolations of the 
curvature coefficients C\ and c 2 as functions of T c /y/ai. Setting T c /^fd\ equal to its 
physical value, cf. m, yields the “physical” set of curvature coefficients 

ci = 0.41 c 2 = 0.53 . (25) 

Since the ratio cr 2 /cri satisfies (tTKl) for all parameter sets used in the interpolation, it is 
expected to satisfy that condition also at the physical point (1251) . On the other hand, 
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this quantity can also be evaluated directly at the physical point, providing a cross¬ 
check of the interpolation procedure. Measuring <7 2 /<7i directly for the parameters CHD 
yields the value cr 2 /cri = 1.37, deviating only slightly from (ffiflh This is adequate for 
the present simplified model, which is only constructed here for comparison purposes. 
Indeed, in the comprehensive model discussed further below, the analogous cross-check 
of the interpolation procedure yields complete agreement with m at the physical point. 

Having determined the “most physical” set of curvature coefficients (GSD within this 
restricted model framework, one can predict the behavior of the string tensions at finite 
temperatures. The most significant quantitative prediction which is accessible in this 
way is the behavior of the spatial quark and diquark string tensions, af and erf, in 
the deconhned phase. For the parameter set (125)) . the deconhned phase is realized for 
N t — 1, which corresponds to the temperature 9 T = 1.9 T c , as one can infer by measuring 
the zero-temperature quark string tension cq in lattice units and using m- At this 
temperature, one obtains 

of(T= 1.9 T c )/cr 1 (T = 0) = 1.61 a§(T = 1.9 T c )/a 2 (T = 0) = 1.46 . (26) 

Thus, the spatial diquark string tension rises less strongly in the deconhned phase than 
the spatial quark string tension; the ratio between the two at T = 1.9 T c is reduced to 

of (T = 1.9T c )/crf(T = 1.9 T c ) = 1.25 (27) 

compared to the zero-temperature value 02/01 = 1.36. 

5.2 Effect of world-surface area term 

As a hrst attempt to generalize the action in order to obtain a more realistic phe¬ 
nomenology, the effect of the world-surface area term m was briehy explored. This 
type of term was already investigated within the SU(2 ) and 577(3) models [21,24] and 
ultimately discarded, since it did not enhance the phenomenological hexibility of those 
models. In the 577(4) model, it can be used to differentiate between the two types of 
vortices without abandoning C\ — c 2 . Thus, the condition C\ = c 2 was maintained, and 
instead m with A(0) = A(l) = 0 but nonvanishing A(2 ) was introduced in the hope 
of maintaining a discernible first-order deconfinement transition while approaching the 
condition (CHI). 

As an example, using A{2) = 0.130, the set of curvature coefficients c\ = c 2 = 0.460 
yields <r 2 /cq = 1.17 while simultaneously realizing the deconfinement transition on a 
N t = 2 lattice, cf. Fig. d The behavior is similar to the one observed above in the 
case of a pure curvature action; approaching condition CHI) weakens the action density 
discontinuity at the phase transition. In fact, in the present case, the discontinuity is 

9 Note that string tensions at other temperatures above T c are not directly accessible at (12511 . but 
can be obtained using interpolations of measurements at other parameter sets to the physical point 
(pf cf. [21]. This was not pursued further here. 
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Figure 4: Distribution of the action density s at the deconfining phase transition, for 
an action containing a vortex world-surface curvature term with c\ = C 2 = 0.460 and a 
world-surface area term with A(2 ) = 0.130. The measurement was taken on a 30 3 x 2 
lattice. 

even more strongly suppressed than in the analogous pure curvature case at the same 
ratio cr 2 /cri = 1.17. Thus, introducing the world-surface area term exacerbates the 
problem rather than alleviating it. 

Increasing A{2) further, when (fTHl) is fully satisfied, the discontinuity indeed has 
disappeared again. When A(2) = 0.300, the deconfinement transition occurs at c\ = 
c 2 = 0.439, and (TTHl) simultaneously holds. The author again carried out a careful scan 
of Ci = c 2 in steps of 0.0001 around this point and found no evidence of first order 
behavior. 

5.3 Branching term 

As has become clear from the preceding sections, a vortex model action based purely 
on the world-surface characteristics curvature and area is no longer phenomenologically 
viable in the SU (4) case. This provides the motivation for introducing the new branch¬ 
ing term m into the vortex model action. The reason for this particular choice lies 
in the following observation. Both the SU(2) and the 577(3) vortex models [21,24] are 
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governed by the same type of action; nevertheless, the 577(3) model exhibits a clear 
first-order deconfinement transition, whereas the SU(2 ) model does not. The difference 
between the two models lies in the class of allowed vortex topologies; only the 57/(3) 
model provides for vortex branching. In view of this, it seems plausible to assume that 
facilitating vortex branching is conducive to a first-order deconfinement phase transi¬ 
tion. Therefore, introducing a new action term encouraging vortex branching, such as 
0, into the SU (4) vortex model is expected to provide a viable mechanism for restor¬ 
ing phenomenologically correct behavior at the phase transition. This expectation is 
confirmed by the measurements discussed in the following. 

The physical point of the model was determined in analogy to the procedure used 
in section 15 .11 Properties of the deconfinement transition generally are not directly 
accessible for the physical set of coupling constants; instead, they are inferred by inter¬ 
polation of phase transition data obtained at different N t . Specifically, for N t = 1, 2, 3 
the parameter sets were determined which realize the deconfinement temperature while 
simultaneously satisfying both of the conditions 0 and (1241) . These parameter sets 
are listed in Table [21 and corresponding action density distributions illustrating the 
first-order behavior of the phase transitions are displayed in Figs. IH1T71 



Cl 

c 2 

b 

Tc/y/(TT 

N t = 1 

0.2785 

0.4005 

0.1403 

0.90 

N t = 2 

0.4558 

0.7983 

0.6950 

0.61 

N t = 3 

0.5925 

0.7059 

0.3800 

0.50 


Table 2: Sets of coupling constants realizing the deconfinement temperature as well as 
satisfying two of the three conditions defining the physical point, namely ED and (121 . 
The remaining condition ED is subsequently satisfied by interpolating the parameters 
using the data in the final column. 

The table furthermore records the ratio T c / yTrj" corresponding to each parameter 
set. To obtain the physical point, one constructs the quadratic interpolations of the 
coupling constants ci, C 2 and b as functions of T c / Setting T c /y/d\ equal to its 
physical value, cf. ED, yields the physical point 

ci = 0.45 c 2 = 0.80 b = 0.71 (28) 

As is evident from Tabled the physical point is very near the parameters obtained for 
N t — 2; the uncertainties inherent in the interpolation thus remain small. Since the 
conditions 0 and ( 121 ) were satisfied for all parameter sets used in the interpolation, 
they are expected to be satisfied at the the physical point 0 as well. In the case of 
0, this can be cross-checked. Measuring <t 2 /<7i directly at the physical point indeed 
yields agreement with 0. Of course, (121 cannot be cross-checked directly at the 
physical point. 
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Figure 5: Distribution of the action density s at the deconfining phase transition, for an 
action containing a vortex world-surface curvature term with c\ = 0.2785, c 2 = 0.4005 
and a branching term with b = 0.1403. The measurement was taken on a 30 3 x 1 lattice. 
Taking the distance between the two maxima as a rough measure of the action density 
discontinuity As and supplementing this with a measurement of the zero-temperature 
quark string tension oy in lattice units at the same set of coupling parameters, the 
condition (1^11) is satisfied by this first-order transition. 


To appreciate the ultraviolet cutoff scale, it is useful to cast the lattice spacing a in 
physical units. Setting the zero-temperature quark string tension to ay = (440 MeV) 2 , 
and combining this with the measurement of that string tension in lattice units at the 
physical point, oya 2 = 0.68, yields 


a = 0.37 fin . (29) 

Within the framework of the random vortex world-surface model, this simultaneously 
characterizes the transverse thickness of the vortices, since it implements a minimal 
distance which parallel vortices must be apart in order to remain distinguishable. In 
the present 577(4) case, this distance is slightly smaller than in the SU(2) and SU(3) 
cases, in which the minimal distance is 0.39 fm. As already noted further above, the fact 
that both of the two physically distinct types of vortices present in the /SC/(4) model are 
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Figure 6: Distribution of the action density s at the deconfining phase transition, for an 
action containing a vortex world-surface curvature term with c\ = 0.4558, C 2 = 0.7983 
and a branching term with b = 0.6950. The measurement was taken on a 18 3 x 2 lattice. 
Taking the distance between the two maxima as a rough measure of the action density 
discontinuity As and supplementing this with a measurement of the zero-temperature 
quark string tension cri in lattice units at the same set of coupling parameters, the 
condition (E£Hl is satisfied by this first-order transition. 


treated as having the same thickness is a model restriction due to the present hypercubic 
lattice formulation; in general, there is no physical reason forcing the two thicknesses to 
be the same. Note also that a decrease of the average vortex thickness as the number 
of colors is raised seems plausible in view of the qualitative picture discussed in section 
12.41 cf. Fig. d 

Having modeled the lattice Yang-Mills data presented in section d it is possible to 
predict string tensions at finite temperatures. In particular, the quantitative behavior 
of the spatial string tensions erf and erf in the deconhned phase can be accessed. At the 
physical point (BED, a lattice with N t = 1 realizes the temperature 10 T = 1.95 T c , and 

10 This temperature determination is minimally adjusted compared to the value quoted in the prelim¬ 
inary report [34] due to improved statistics in the calculation of the zero-temperature string tension. 
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Figure 7: Distribution of the action density s at the deconfining phase transition, for an 
action containing a vortex world-surface curvature term with ci = 0.5925, C 2 = 0.7059 
and a branching term with b = 0.3800. The measurement was taken on a 20 3 x 3 lattice. 
Taking the distance between the two maxima as a rough measure of the action density 
discontinuity As and supplementing this with a measurement of the zero-temperature 
quark string tension cri in lattice units at the same set of coupling parameters, the 
condition (E£Hl is satisfied by this first-order transition. 

one obtains 

crf(T = 1.95 T c )/cr 1 (T = 0) = 1.34 af(T = 1.95 T c )/a 2 {T = 0) = 1.44 . (30) 

While both string tensions exhibit the characteristic rise with temperature in the decon- 
fined phase, the behavior of erf in particular is quite different from the simplified model 
discussed in section o In contradistinction to that model, here, erf rises less strongly 
than erf, and the ratio of the spatial diquark string tension to the spatial quark string 
tension is enhanced at T = 1.95 T c , 

of (T = 1.95 T c )/erf (T = 1.95 T c ) = 1.46 (31) 

compared to the zero-temperature value cr 2 /cr x = 1.36. This should be contrasted with 
the behavior in the aforementioned simplified model, where of(T = l-9T c )/of (T = 
1.9 T c ) = 1.25, cf. fl23). 


20 



Furthermore, in the more recent survey of SU( 4) lattice Yang-Mills characteristics 
[33], which became available to the author after the bulk of the numerical work in the 
present investigation was completed, data for this string tension ratio are reported. 
According to these data, the ratio erf /erf remains near its zero-temperature value in 
the temperature regime considered here. Thus, while the comprehensive vortex model 
constructed in the present section adjusts the value of erf/erf in the correct direction 
compared to the simplified model of section 15.11 the adjustment overshoots the SU(4) 
lattice Yang-Mills result considerably. The result m is almost as far above the lattice 
data as the corresponding result in the simplified model, ( 22 ), is below those data. 

ft should be noted that readjusting the present model to reproduce updated values 
from [33] for the observables discussed in section Q] would of course impact this com¬ 
parison. On the other hand, regardless of the outcome of such an exercise, a further 
refinement of the SU(4) vortex model action would presumably permit fitting also the 
finite-temperature spatial string tension ratio correctly, in addition to the observables of 
section El For instance, one could contemplate using a coupling constant c 12 ^ Ci • c 2 in 
the curvature action (USD; while this option was discarded above due to its physical sim¬ 
ilarity to the branching action m, it is not precisely the same and thus may allow for 
some additional quantitative tuning. Another possible source of discrepancies between 
the model constructed here and SU (4) Yang-Mills theory is the artificial restriction to a 
common thickness for both types of vortices present for SU (4) color. Refinements along 
these lines will not be considered further here. If anything, such an endeavor would only 
further reinforce the conclusions drawn in the next, concluding, section below. 


6 Conclusions 

The goal of this investigation was to construct a random vortex world-surface model 
for S77(4) Yang-Mills theory, particularly with a view towards the question whether 
infrared S77(4) Yang-Mills phenomenology forces one to abandon the dynamical concept 
employed in the SU(2 ) and SU( 3) cases, which is based purely on vortex world-surface 
characteristics. General arguments, put forward in [30] and reviewed in section|21 suggest 
that a shift away from this type of dynamics occurs as the number of colors is raised, 
with additional color structures present on the vortex world-surfaces attaining their own 
dynamical significance and influencing the vortex ensemble. 

The modeling effort carried out in this work supports these arguments. It proved im¬ 
possible to construct a vortex model which reproduces the main infrared characteristics 
of SU(4) Yang-Mills theory based purely on a vortex world-surface curvature action. 
Instead, an explicit weighting of vortex branchings, which, depending on the color de¬ 
scription, are associated with magnetic monopoles or nexi, turned out to furnish the 
additional flexibility needed to match SU( 4) Yang-Mills theory. 

The need to introduce new action terms of course impacts the predictive capability 
of the model. Indeed, the effortless predictivity of the SU(2 ) and 577(3) models, in 
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which merely one dimensionless curvature coefficient needed to be adjusted, is lost in 
the 577(4) case. Three dimensionless parameters were necessary to obtain a sufficiently 
flexible model. After having determined these parameters, the behavior of the spatial 
string tensions at high temperatures was predicted and compared to newer lattice Yang- 
Mills data which became available subsequently, ft should be noted that the predictions 
are obtained at the ultraviolet limit of validity of the vortex model, and thus can be ex¬ 
pected to be more prone to error than observables at lower temperatures; nevertheless, 
the discrepancy in the ratio of the spatial diquark to the spatial quark string tension 
at T « 2T c , amounting to about 7%, seems significant. Aside from the possibility that 
this comparison may improve after readjusting the present vortex model to reproduce 
updated values from [33] for the observables discussed in section El further model re¬ 
finements, as suggested at the end of section 15.31 would presumably permit closing this 
gap in any case, ffowever, this would merely exacerbate the loss of predictive capability, 
already noted above, compared with the SU{2) and S77(3) vortex models. 

In view of this, it does not seem attractive to pursue the construction of a yet 
more detailed S77(4) vortex model. Indeed, the main question which this investigation 
aimed to address can be answered in the affirmative already at the present stage, as 
discussed above: As the number of colors is raised beyond N — 3 in SU(N) Yang-Mills 
theory, the corresponding infrared effective vortex description cannot continue to rely on 
dynamics based purely on vortex world-surface characteristics. Solid evidence for this 
emerges in the N = 4 case studied in this work. To restore the correct 577 (4) Yang-Mills 
phenomenology at the deconfining phase transition, a generalization of the vortex model 
action was necessary, and an explicit dependence of the dynamics on vortex branching 
characteristics proved sufficient to achieve that goal. 
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